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Abstract 

The spin-lattice relaxation rate 1/Ti and the spin echo decay rate \/T2g 
for the spin-^ antiferromagnetic Heisenberg chain are calculated using quan- 
tum Monte Carlo and maximum entropy analytic continuation. The results 
are compared with recent analytical calculations by Sachdev. If the nuclear 
hyperfine form factor Aq is strongly peaked around q = ir the predicted low- 
temperature behavior [1/Ti ~ ln^^^{l/T), I/T2G ~ In'^/^il/T) /VT] extends 
up to temperatures as high as T/J ~ 0.5. If Aq has significant weight for 
(7 ~ there are large contributions from diffusive long-wavelength processes 
not taken into account in the theory, and very low temperatures are needed 
in order to observe the asymptotic T — > forms. 
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The one-dimensional spin-| antiferromagnetic Heisenberg hamiltonian, 

H = JJ2S.,-S.,+,, (1) 

i 

is relevant as a starting point for understanding the magnetic properties of many quasi 
one-dimensional materials. Examples include CuCl-2N(C5H5),0 KCuFsP and several tetra- 
cyanoquinodimethan (TCNQ) charge transfer salts-iS NMR and NQR are commonly used 
techniques for studying the spin dynamics of materials such as those listed above. The low- 
frequency dynamic susceptibility is accessible through the spin-lattice relaxation rate 1/Ti 
and the spin echo decay rate I/T2G. Theoretical results for the temperature dependence 
of both these rates were recently obtained by Sachdev,!! using a form for the dynamic sus- 
ceptibility first derived by Schulzi using the Bosonization method. Neglecting logarithmic 
corrections, 1/Ti is predicted to be constant at low temperature, and I/T2G is predicted 
to diverge as T~^/^. With logarithmic corrections taken into account both rates acquire a 
factor ln^/^(l/T). These results are expected to be valid only for temperatures T <^ J, and 
it is important to verify their validity as well as to obtain results also at higher tempera- 
tures. Here results are presented for 1/Ti and I/T2G computed using quantum Monte Carlo 
(QMC) simulations of chains of up to 1024 spins. I/T2G is related to static susceptibilities 
directly computable in the simulations. The dynamic susceptibility required for extracting 
1/Ti is calculated in imaginary time and continued to real frequency using the maximum 
entropy method.0i 

The results for the temperature dependence of both 1/Ti and I/T2G at low temperatures 
are in good agreement with Sachdev's predictions. At higher temperatures diffusive modes 
not taken into account in the theory cause significant deviations. If the nuclear hyperfine 
form factor has large weight at long wavelengths very low temperatures are needed for the 
asymptotic forms to apply, and they may then be difficult to observe experimentally. The 
results presented here should be useful for comparisons with experiments also at higher 
temperatures. 

The NMR spin-lattice relaxation rate is given byi 



;^ = iEEKI'^(5>^-o), (2) 

-'l a q 

where A'^ is the hyperfine form factor, and a denotes the two axes perpendicular to the 
external field direction. S'(g, u) is the dynamic structure factor, which is related to the 
imaginary part of the dynamic spin susceptibility according to S{q, uj) = x"{q^ ^) I (1 — e~^'^) , 
where /3 = l/ksT. Here an isotropic form factor = Aq will be assumed. Defining 

^^M = iEiAr^(?'^)' (3) 
<? 

the spin-lattice relaxation rate is then obtained as 1/Ti = 2SA{i^ — * 0). For the numeri- 
cal calculations carried out here it is more convenient to work directly with the hyperfine 
coupling A{r) in coordinate space. Define 

Ca{t) = J2A{r,)A{r,){Sar)S^m, (4) 

where S'f (r) = e'^^ Sfe~'^^ . With C^(t) calculated numerically, Sa{^^) can be obtained by 
inverting the relation 

oo 

C^(r) = - / d^5A(a;)e-™ (5) 

TT J 
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using the maximum entropy technique.0i This method is described in detail in Ref. |^, and 
was recently applied in a calculation of the spin-lattice relaxation rate of the two-dimensional 
Heisenberg model. 

The gaussian component of the spin echo decay rate is related to the the nuclear spin- 
spin interactions mediated by the electrons. Under conditions discussed by Pennington and 
Shchtei0 

7^=[lT.J'iO,xf\ (6) 
where Jzixi, X2) is the ^-component of the induced interaction between nuclei at xi and X2: 
Jzixi, X2) = ^(^1 - ri)A{x2 - rj)x{i - j). (7) 



The static susceptibility ~ j) is given by the Kubo formula 

x{t'3) = ldT{s:iT)s]m. (8) 



The hyperfine interaction A{r) is normally very short ranged. Here a situation is con- 
sidered where the nuclei studied reside at the sites of the electronic spins modeled by the 
hamiltonian (P. The hyperfine coupling is assumed to have a direct contact term of strength 
A{0), and a transferred nearest-neighbor term of strength A{1). Results are presented for 
several values of the ratio R = A{1)/A{0). 

A stable inversion of the relation requires that Ca{t) is known to very high accuracy. 
Here a quantum Monte Carlo method based on stochastic series expansionS (a generaliza- 
tion of Handscomb's methodEl) is used. This technique is free from systematical errors of the 
"Trotter break-up" used in standard methods.0 The imaginary time correlation functions 
needed have been calculated to within relative statistical errors of only 10~^ or lower for tem- 
peratures down to T = J/8. This high accuracy is required for obtaining a reliable estimate 
of Sa{uj 0). The static susceptibilities are computed directly in the QMC simulation, 
and hence the calculation of I/T2G is not hampered by potential problems associated with 
analytic continuation. Accurate results for I/T2G have been obtained at temperatures as 
low as T = J/32 for systems of up to 1024 sites, which is large enough for finite-size effects 
to be completely negligible. 

In order to test the accuracy of a calculation of 1/Ti by analytic continuation of QMC 
data, complete diagonalizations of 16-site chains were also carried out. Comparisons of 
S{uj) = jj: J2q S{q, uj) obtained in these calculations with numerically continued QMC data 
are shown in Fig. 1. The maximum entropy method requires a "default model" which defines 
the zero of entropy.Bi In all calculations presented here a flat default model was used. Exact 
diagonalization gives >S'(ti;) as a finite number of delta functions and their corresponding 
weights. Here the results are plotted as histograms in order to facilitate comparison with 
the maximum entropy result. The jagged structure of the exact diagonalization result, 
which is due to the small size of the system, cannot be reproduced by the maximum entropy 
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method. The results do however represent reasonable frequency averages. Note that even the 
high-frequency behavior is obtained quite accurately. Clearly a 16-site system is not large 
enough for extracting the low- frequency behavior at low temperatures. The 1/Ti results 
presented below are for systems of 256 spins, and the accuracy of the imaginary time data 
used for the analytic continuation is even higher than the data used for the 16-site results 
shown in Fig. 1. Comparisons with results obtained for 128 spins indicate that there are no 
significant finite size effects at the temperatures considered. 

Both fluctuations of the uniform (g ^ 0) and staggered (g ~ tt) magnetization con- 
tribute to the NMR rates of half-integer spin chains. At low temperatures the staggered 
contribution dominates.!! Neglecting the uniform fluctuations, Sachdev obtained the asymp- 
totic low-temperature forms (in units where ^ = /c^ = 1) 

Ti c 

1 r~r 

where / ~ 8.4425, and c is the spinon velocity, which for spin-i is c = |. D is the prefactor 
of the asymptotic equal-time spin-correlation function, which is not known accurately. 0The 
marginally irrelevant operator present for the critical spin chains has not been taken into 
account in the derivation of the above forms. This is expected to lead to a multiplicative 
correction ln^/^(A/T) for both 1/Ti and l/T2G'.i Hence, the ratio T2g/{VTTi) should be a 
constant, even with logarithmic corrections included. 

If I/T2G ~ ln^/^(A/T)/-\/T as predicted by Sachdev, T(1/T2g)^ should be a linear func- 
tion of ln( J/T). In Fig. 2, T{AIT2g)~^ is graphed versus In (J/T) for several values of the 
hyperfine coupling ratio R = A{1)/A{0). In cases where the corresponding Aq is peaked 
around q = n {R < 0) a linear behavior is seen in a wide temperature regime. The points 
for R = —0.25 and R = —0.5 nearly coincide at low T, indicating that the q ^ tt contribu- 
tions almost completely dominate the behavior in both cases. A line fit to the R = —0.5 
points gives A = 0.92 J and the amplitude D = 0.080 in (^^. For R > contributions 
from q ^ rapidly become important at high temperatures, and for large values of R the 



asymptotic behavior can only be observed at very low temperatures. 

Results for 1/Ti divided by are shown in Fig. 3. The expected weak (logarithmic) 
increase as T decreases can be seen below T/J = 0.5 if R is large and negative, so that 
q ^ TT processes dominate Sa{uj — > 0) even at relatively high temperatures. For R> 0, 1/Ti 
decreases with decreasing T down to quite low temperatures — for R = 0.25 this behavior 
extends down to the lowest temperature studied. The enhancement of 1/Ti at high T is 
caused by the diffusive g ~ processes not taken into account in the forms In order to 
more clearly determine the importance of these modes one can study the ratio 

Sq<n/2{^ 0) _ T,q<7T/2S{q,UJ 0) 

which is graphed versus the temperature in Fig. 4 (these calculations were carried out on 
systems of 128 spins). At T = J the q < tt/2 contribution is approximately 50%, and 
decreases rapidly at lower temperatures. These results confirm Sachdev's conclusion! that 
the g ~ contribution to 1/Ti is negligible in the limit T 0. 

Returning now to the results shown in Fig. 3, there are not enough low-temperature 
data to extract the asymptotic temperature dependence of 1/Ti. The results are, however, 
consistent with a divergence of the predicted form ln^/^(A/T) with the same A = 0.92J as 
was found above for I/T2G. The amplitude needed in Eq. (^bj) is then D ^ 0.14, which 
is significantly larger than the amplitude extracted from I/T2G above. Hence, the ratio 
T2g/{^/TTi) is different from Sachdev's prediction. The ratio is graphed versus temperature 
in Fig. 5. For i? < it is indeed almost constant below T / J 0.5, whereas for i? > there 
is a significant temperature dependence down to the lowest temperatures considered. For 
positive R there is a sharp maximum in T2G / {'s/TTi) , arising from the minimum in I/T2G 
seen in Fig. 2. The R = —0.5 result for T2g/{'\/TTi) at low temperatures is approximately 
3.0 — 3.1, which is almost a factor 2 larger than what is obtained from Eqs. (^. 

In summary, the NMR rates 1/Ti and I/T2G have been calculated for the spin-1/2 
Heisenberg model, using quantum Monte Carlo and maximum entropy analytic continua- 
tion. The temperature dependence at low temperature is in good agreement with Sachdev's 
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recent theoretical results, which include only the contributions from staggered magneti- 
zation fluctuations. At high temperature damped q ^ modes are important, and can 
dominate the NMR rates if the hyperfine form factor has large weight at long wavelengths. 
In such cases very low temperatures are needed to observe the asymptotic forms. In many 
real systems effects of interchain couplings may become important before the asymptotic 
regime is reached, and the low-temperature forms may therefore not be easily observed. 
The results here should then be useful for determining the relevance of a description by the 
one-dimensional Heisenberg model based upon measurements at higher temperatures. It 
can be noted that early NMR experimentsi on (NMP)(TCNQ) indicate a behavior of 1/Ti 
similar to the result shown in Fig. 3 for a small positive hyperfine ratio A{1)/A{0), with no 
indication of a low-temperature increase down to T ^ O.IJ. 

It will be interesting to apply the techniques used here to calculate the NMR rates of other 
one-dimensional systems. Work on coupled spin chains is in progress.lll Itinerant electrons 
described by one-dimensional Hubbard-type models, including electron-phonon interactions, 
can also be studied. 
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FIGURES 

FIG. 1. QMC and exact diagonalization results for S{u;) — j^J2gS{q,u;) of a 16-site 
chain. The exact results are plotted as histograms, and the curves are numerically continued 
QMC imaginary time correlation functions. 

FIG. 2. Results for the spin echo decay rate, graphed as T{AIT2g)^'^ vs. ln( J/T) for 
different hyperfine ratios R = A{1)/A{0). 

FIG. 3. Spin-lattice relaxation rates vs. temperature for different hyperfine ratios R — 
A{1)/A{0). The solid curve is of the form ln(A/T), with A = 0.92 J. 

FIG. 4. The long- wavelength contribution to — > 0) = J2qS{q,u; — > 0) vs. tempera- 
ture. 

FIG. 5. The ratio T2g/VTTi vs. T for different hyperfine ratios R = A{1)/A{0). 



9 




CO 



FIG. 1, A. W. Sandvik 




ln(J/T) 



FIG. 2, A. W. Sandvik 



A 



o R=0 
■ R=-0.25 
□ R=-0.50 
A R=+0.10 
A R=+0.25 




o o 



□ □ □ □ 



0.0 



0.2 



0.4 



0.6 



T/J 



0.8 



1.0 
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